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I.  INTRODUCTION 


Interior  ballistics  flows  of  conventional  charges  in  gun  tubes  possess 
complex  flow  patterns.  Their  complexity  is  due  to  both  the  heterogeneous 
structure  of  the  charge  and  the  rapidly  changing  flow  conditions  within  a  few 
milliseconds.  The  fast  rise  in  pressure  and  temperature  caused  by  the  burning 
of  the  propellant  initiates  a  turbulent,  multidimensional,  multiphase  flow 
which  is  coupled  with  the  accelerating  projectile  motion.  A  complete 
mathematical  model  that  describes  all  the  physical  phenomena  occurring  in  an 
interior  ballistics  cycle  is  not  presently  available.  However,  several  models 
that  simulate  some  of  the  phenomena  exist  or  are  being  developed. * 

On  the  other  hand,  it  is  not  possible  for  technical  reasons  to  make 
detailed  experimental  measurements  of  the  complete  interior  ballistic  cycle. 
Some  standard  techniques  as  well  as  some  new  special  techniques  under 
development  determine  only  specific  quantities  in  real  weapons  or  in 
simulators  under  simplified  flow  conditions.  Commonly  measured  quantities  are 
the  gas  pressure  and  projectile  motion.  Other  quantities  such  as  the 
temperature  distribution  in  the  gas  and  in  the  gun  tube  wall,  the  velocity 
distribution  of  the  gas  and  solid  particles  inside  or  outside  of  boundary 
layers,  the  particle  distribution,  the  turbulence  pattern,  etc.,  cannot  be 
accurately  determined  by  experiment.  Thus,  a  need  exists  for  modelling  of  the 
interior  ballistic  cycle  so  that  the  dynamic  development  of  these  quantities 
can  be  studied,  and  their  impact  on  ballistic  problems  can  be  evaluated. 

A  new  computational  capability  for  the  investigation  of  interior 
ballistics  flows  is  the  DELTA  code,  which  is  under  development  at  the 
Ballistic  Research  laboratory.  The  purpose  of  this  code  is  to  address 
particular  ballistics  problems  related  to  the  boundary  layer  development,  the 
heat  transfer  to  the  tube  wall,  the  "urbulence,  and  the  time-dependent 
distribution  of  additive  particles.  In  the  following,  we  shall  give  a  short 
description  of  the  DELTA  code  and  present  computational  results  of  an 
investigation  of  heat  conduction  and  turbulence  effects. 


II.  REVIEW  OF  THE  MODELS  IN  DELTA 


The  flow  which  is  modeled  by  DELTA  is  a  multidimensional,  two-phase  flow 
inside  a  gun  tube.  Presently,  the  flow  is  assumed  to  be  axisymmet r ic .  At  the 
rear  end,  the  so-called  breech,  the  tube  is  closed  by  a  stationary  flat  plate 
while  the  front  boundary  is  a  moving  flat-based  projectile.  The  flow  is 
assumed  to  be  viscous  and  heat  conducting,  and  it  can  be  either  laminar  or 


^ "Fluid  Dynamics  Aspects  of  Internal  Ballistics,"  AGARD  Advisory  Report 
No.  172,  1982. 


turbulent.  The  wall  may  be  either  adiabatic  or  it  may  allow  heat  transfer 
from  the  gas.  Heat  transfer  is  restricted  to  the  tube  wall  (excluding  the 
breech).  The  core  flow  is  fully  coupled  to  the  moving  projectile,  to  the 
boundary  layer  development  and,  if  desired,  to  the  heat  conduction  in  the  tube 
wall.  By  fully  coupled  we  mean  that  each  of  these  phenomena  can  affect  all 
others.  For  example,  the  boundary  layer  development  can  alter  the  details  of 
the  core  flow.  This  would  not  be  the  case,  for  example,  in  a  boundary  layer 
type  model. 

The  mathematical  model  in  DELTA  (the  balance  equations  for  the  gas-phase 
and  one  solid  phase),  is  based  on  an  unsteady  volume-averaged  formulation. 

The  gas  phase  is  described  by  averaged  equations  corresponding  to  the  full 
Navier-Stokes  equations  for  a  compressible  fluid.  The  model  is  closed  by 
averaged  coefficients  of  viscosity  and  heat  conduction,  an  averaged  viscous 
stress  tensor,  an  averaged  dissipation  function  and  an  averaged  heat 
conduction  function.  Since  interior  ballistics  flows  usually  produce  high  gas 
pressures,  the  Noble-Abel  equation  of  state  is  used  so  that  some  real  gas 
effects  can  be  included.  The  gas  turbulence  is  represented  by  algebraic 
mixing  length  models.  The  solid  phase  is  described  by  the  averaged  equations 
for  arrays  of  incompressible  particles  which  can  undergo  deformations.  The 
derivation  of  the  equations  is  given  in  Ref.  2.  The  equations  of  the  solution 
algorithms  are  listed  in  Ref.  3. 

The  system  of  partial  differential  equations  for  the  ax i symmetric  two- 
phase  flow  region  is  solved  by  a  linearized  Alternating  Direction  Implicit 
(ADI)  scheme.  This  scheme  transforms  the  differential  equations  into  a  system 
of  linear  algebraic  equations.  The  corresponding  matrix  has  a  block, 
tridiagonal  structure,  allowing  an  efficient  determination  of  the  solution  at 
each  new  time-step.  Details  about  the  derivation  of  the  scheme  are  presented 
in  Ref.  3. 


III.  HEAT  TRANSFER  TO  AND  TEMPERATURE  DISTRIBUTION  IN  THE  TUBE  WALL 

The  heating  of  the  gun  tube  wall  caused  by  convection,  heat  conduction, 
and  radiation  of  the  hot  propellant  gas  enhances  the  gun  tube  wear  and 
erosion,  and  therefore,  affects  the  lifetime  of  gun  tubes.  An  experimental 
determination  of  the  inner  wall  surface  temperature  is  quite  difficult. 
Commonly  used  thermocouples  ire  limited  use  in  interior  ballistic 


2 

“f. eic.it. s,  A.K.F.,  Sc.nmitt,  d.A..  "To ree-Dimensiona  1  Model in*  of  as- 
Combusting  Solid  Two-Phase  Flows,  "Multi-Phase  Flow  and  Heat  Transfer 
III,  Part  B:  Applications,  T.  N,  Veziroglu  and  A.E.  Bergles,  editors, 
pp.  581-598,  Elsevier  Science  Publishers,  Amsterdam,  1984. 

jSchiritt,  d.A.,  "A  .utre  ri  ca  t  Algorithm  'or  the  Multidimensional,  Multiohase. 
Viscous  Equations  of  Interior  Ballistics,"  Transactions  of  the  Second  Army 
Conference  on  Applied  Mathematics  and  Computing,  ARO-Report  85-1, 
pp.  849-59 1 ,  1985 . 


applications  because  such  applications  require  accurate  measurements  in  very 
short  time  intervals,  a  close  contact  to  the  flow,  and  special  thermal 
properties  of  the  gauges. 


There  are  several  mathematical  models  of  the  heat  transfer  to  the  gun 
tube  wall.  They  can  be  divided  in  four  categories  according  to  their 
complexity.  The  first  type  uses  a  very  simplistic  boundary  layer  calculation 
and  a  hea£  transfer  correlation,  e.g.,  Colburn's  analogy,  to  obtain  the  heat 
transfer.14  Heat  losses  in  the  core  flow  are  considered.  The  empirical  heat 
transfer  correlations  are  derived  for  fully  developed,  steady,  one-phase  pipe 
flow.  The  main  feature  of  this  type  of  model  is  the  emphasis  on  the 
calculation  of  the  core  flow.  In  models  of  the  second  category  the  emphasis 
is  on  a  description  of  the  boundary  layer  by  using  more  general  boundary  layer 
equations.  "  The  heat  transfer  to  the  tube  wall  again  is  described  by 
correlations.  The  boundary  layer  edge  is  not  coupled  with  the  computation  of 
the  core  flow.  Instead,  one  assumes  for  the  latter  values  which  represent 
approximately  the  core  flow.  The  third  category  makes  use  of  general  boundary 
layer  equations  and  of  a  balance  of  heat  fluxes  from  the  hot  gas  to  the  gun 
tube  wall  at  the  inner  wall  surface. ^  The  conditions  at  the  boundary  layer 
edge  are  comparable  to  those  in  the  second  category.  Neither  of  these 
approaches  includes  all  the  feedback  mechanisms  from  the  tube  to  the  core 
flow,  and  therefore,  to  the  projectile  motion.  The  fourth  category  of  models 
uses  a  fully  coupled  approach  whereby  the  phenomena  in  the  core  flow  are 
directly  linked  to  the  projectile  motion,  the  boundary  layer  development,  the 
heat  transfer  to  and  the  heat  conduction  in  the  tube  wall,  and  vice  versa. 

This  is  achieved  by  using  a  single  system  of  equations  everywhere  in  the  gas 
region  and  a  coupled  system  of  heat  conduction  equations  in  the  tube  wall. 

The  solution  of  these  sets  automatically  provides  the  boundary  layer  solution 
in  the  boundary  layer  region,  the  core  flow  solution  in  the  core  flow,  and  all 
the  necessary  coupling  that  naturally  occurs  In  the  flow.  Although  this  type 


^Shelton,  S. ,  Bergles,  A.,  Saha,  ?.,  "Study  of  Heat  Transfer  and  Erosion  in 
Gun  Barrels,"  Air  Force  Armament  Laboratory,  Eglin  Air  Force  Base,  Florida, 

AF  AT  L-TR-7  3-69 ,  1973 . 

'’Gough,  P. ,  "Modeling  of  Rigidized  Gun  Propelling  Charges,"  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  Maryland  ARBRL-CR-00318,  J  9  8  3 . 
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WaLL,"  J.  BaLlistics  6  (3),  pp.  1518-1524,  1982. 

7£arlett,  E.P.,  Anderson,  L..W. ,  Kendall,  R.H.,  'Time-Dependent  Boundary  Layers 
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of  solution  is  most  complex,  it  provides  a  flow  description  with  the  fewest 
assumptions  and  approximations.  In  light  of  the  scarcity  of  experimental 
measurements  with  which  to  compare  the  calculations,  we  feel  such  an  approach 
to  be  best.  The  DELTA  code  is  an  example  of  a  fourth  category  model. 


The  heat  transfer  model  in  DELTA  consists  of  the  equations  governing  the 
heat  conduction  in  the  tube  wall,  and  the  boundary  conditions  which  couple  the 
temperature  in  the  wall  to  the  flow  inside  and  outside  the  tube.  The  heat 
conduction  in  the  tube  wall  is  described  by  the  two-dimensional,  nonlinear 
axisymmetric  equation  for  the  wall  temperature  Tw(t,r,z): 
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The  variables  t,  z  and  r  denote  the  time,  axial  coordinate  in  the  wall  and 
radial  coordinate  in  the  wall,  respectively.  The  specific  heat  cw  and  the 
thermal  conductivity  Aw  of  gun  tube  steel  strongly  depend  on  the 
temperature. 10  By  coraparision,  the  density  of  the  steel  pw  varies  only  a 
little  with  temperature. 

The  most  important  boundary  condition  is  at  the  inner  tube  wall  surface 
where  the  coupling  of  the  flow  region  and  the  wall  occurs.  It  consists  of  the 
balance  of  heat  fluxes  with  a  radiation  effect 


g  3r 


3T  ,  , 

^  +  (.  o  (T4  -TA)  = 


,  5, 

w  3r 


and  of  the  temperature  equilibrium  equation 
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The  variables  A  ,  ,  o ,  T  and  T  denote  the  thermal  conductivity  of 

g  g  g 

the  gas,  the  emissivity  of  the  wall  surface,  the  Stefan-Boltzmann  constant, 
the  gas  temperature,  and  the  maximum  gas  temperature  in  a  given  cross-section 
(z  =  constant),  respectively.  We  emphasize  that  in  DELTA  both  conditions  are 
used  only  at  the  inner  tube  wail  surface,  excluding  the  breech.  The  left  hand 
side  of  the  first  condition  represents  the  heat  flux  on  the  gas  side  towards 
the  wall,  whereby  the  first  term  represents  heat  conduction,  and  the  second 
term  represents  heat  radiation.  The  right  hand  side  gives  Lite  heat  flux  into 
the  tube  wall.  For  the  boundary  conditions  at  the  outer  tube  wall  surface  we 
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chose  the  simple  engineering  condition 
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where  h  is  a  heat  transfer  coefficient  and  Tam^  is  an  outer  ambient 
temperature.  A  more  sophisticated  condition  is  not  needed  (at  least  for  a 
single  shot  weapon)  because  the  heat  usually  does  not  reach  the  outer  surface 
during  a  ballistic  cycle.  Two  additional  boundary  conditions  are  needed  in 
axial  direction.  At  the  projectile  base,  we  set 
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across  the  wall  thickness,  that  is,  we  assume  that  the  projectile  moves  into 
an  area  which  is  at  ambient  temperature.  At  the  breech  and  at  the  projectile 
base,  an  adiabatic  condition 


was  assumed  to  be  adequate. 

The  equations  governing  the  temperature  distribution  in  the  tube  wall  are 
solved  using  the  same  linearized  ADI  method  as  for  the  equations  in  the  gas 
flow  region;  that  is,  the  equations  are  linearized  in  time,  and  are  split 
along  coordinate  directions.  At  each  new  time  level,  we  first  update  the 
temperature  distribution  in  the  wall,  and  then  update  the  dependent  variables 
in  the  flow  region.  This  is  performed  by  the  following  sequence  of  sweeps 
along  coordinate  directions:  an  axial  sweep  followed  by  a  radial  sweep  in  the 
wall,  a  radial  sweep  followed  by  an  axial  sweep  in  the  gas  region,  and  finally 
an  adjustment  of  the  dependent  variables  along  the  inner  wall  surface  to  the 
flux  boundary  condition.  We  omit  a  discussion  of  the  details  of  the  numerical 
procedure  because  they  are  discussed  in  Ref.  3.  In  the  DELTA  code,  the 
thermodynamic  dependent  variables  in  the  gas  region  are  the  specific  gas 
entropy  (s)  and  the  logarithm  of  the  gas  pressure  (q).  Therefore,  the  heat 
flux  boundary  condition  must  be  reformulated  in  terms  of  s  and  q  at  the  new 
unknown  time  level  for  the  radial  sweep  in  the  gas  region.  To  this  end,  we 
transform  the  heat  flux  term  on  the  right  hand  side  of  the  condition  via  the 
chain  rule  as  follows : 
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'•'here  T  is  the  gas  temperature,  and  T  ,  denote  the  partial  deriviatives  of 
f  with  respect  to  s  and  q,  respectively.  H(For  simplicity,  we  dropped  the 
index  g.)  The  linearization  in  time  gives  a  relation  between  the  unknown  new 
time  level  (n)  and  the  known  current  time  level  (c) 
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The  time-deri viati ves  we  approximate  by 
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With  this  approximation,  the  boundary  condition  provides  a  linear  equation  for 
sn  and  qn.  The  equation  is  compatible  with  the  set  of  finite  differenced  and 
linearized  flow  equations. 


IV.  TURBULENCE  MODEL 


In  order  to  estimate  the  influence  of  turbulence  on  the  flow  pattern,  we 
used  two  turbulence  models.  Both  were  equilibrium  algebraic  eddy  viscosity 
models  based  on  Prandtl's  mixing  length  hypothesis.  In  these  models,  a 
turbulent  eddy  viscosity  p  and  a  turbulent  thermal  conductivity  A  are  added 
to  the  molecular  viscosity  p  and  thermal  conductivity  A  ,  respectively, 
yielding  the  effective  values 
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which  then  are  used  in  the  laminar  flow  equations.  The  two  models  differ  in 
the  underlying  assumption  that  the  boundary  layer  consists  either  of  one 
region  or  is  composed  of  two  regions. 


The  one-layer  model  expresses  the  turbulent  eddy  viscosity  by 
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where  P  is  the  local  density,  £  is  Prandtl's  mixing  length,  w  and  u  are  the 
velocity  in  axial  and  radial  direction,  respectively.^  In  the  DELTA  test 
calculations  the  velocity  gradients  were  obtained  from  the  solution  of  the 
governing  differential  equations.  The  mixing  length  was  obtained  from  a 
correlation.  For  a  steady  incompressible  flow  in  a  tube  Nikuradse** 
experimentally  determined  that 

£  =  R  l  0.14  -  0.08  (1  -  £)2  -  0.06  (1  -  £  )4J. 

R  R 


where  R  is  the  tube  radius  and  y  is  the  distance  from  the  tube  wall.  Using 
this  correlation  one  assumes  that  it  models  the  turbulence  also  in  an  unsteady 
and  compressible  flow.  We  test  this  assumption  by  comparing  the  computed 
results  with  those  obtained  using  another  turbulence  model,  a  so-called  two- 
layer  model. 

The  two-layer  model  separates  the  boundary  layer  in  an  inner  and  outer 
region  with  different  formulations  for  each  region.  The  expression 

for  the  eddy  viscosity  in  the  inner  region  is 
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which  is  the  same  as  the  one-layer  formulation.  The  difference  is  in  the 
definition  of  the  mixing  length  £,  which  now  is  calculated  by 


£  =  k  y  D, 

where  k  =  0.4  is  the  von  Karraan  constant,  y  is  the  distance  from  the  wall  and 
D  is  the  van  Driest  damping  factor.  The  latter  is  given  by 
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where 
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and  A  =  26  is  the  van  Driest  constant.  The  subscript  w  indicates  that  the 
subscripted  quantities  are  to  be  evaluated  at  the  wall  surface  (y=0).  The 
wall  shear  stress  x  is  expressed  by 
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In  the  outer  layer  we  use  for  the  eddy  viscosity  the  correlation 
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where  wg  denotes  the  axial  edge  velocity,  6  the  kinematic  boundary  layer 
displacement  thickness  and  6  the  boundary  layer  thickness.  In  our  case  we  is 
the  maximum  axial  velocity  in  the  cross-section  z=constant,  that  is,  the 
velocity  on  the  axis  of  symmetry.  Hence  the  outer  layer  encompasses  in  this 
model  the  core  flow.  The  complete  two-layer  eddy  viscosity  is  given  by 
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turbulent  thermal  conductivity  is  calculated  in  both  cases 
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where  Pr_  =  0.9  is  the  turbulent  Prandtl  number  and  c  is  the  specific  heat  at 
t  P 

constant  pressure. 


V.  RESULTS 


Computational  results  are  presented  for  two  different  types  of  interior 
ballistics  flows.  The  first  type  simulates  a  pure  gas  expansion  flow  behind  a 
projectile  moving  in  a  constant  cross-section  tube.  The  tube  is  closed  at  one 
end  by  a  stationary  surface  cal  Led  the  breech,  and  at  the  other  end  by  the 
movable  flat  based  projectile.  The  initial  states  of  the  gas  are  uniform  and 
quiescent.  Geometrical  data,  initial  conditions  as  well  as  the  thermodynamic 
properties  of  the  gas  are  listed  in  Table  I.  We  designate  this  idealization 
as  the  Lagrange  gun. 

TABLE  I.  LAGRANGE  GUN  PARAMETERS 


Bore  Diameter 

Tube  Length 

Chamber  Length 

Projectile  Mass 

Ratio  of  Specific  Heats 

Covolume 

Molar  Mass 

Initial  Gas  Pressure 
Initial  Gas  Temperature 
Initial  Velocities 


The  flow  in  the  Lagrange  gun  is  very  well  suited  for  the  study  of  several 
important  features  such  as  the  performance  of  the  numerical  procedure,  the 
boundary  layer  development,  the  laminar  and  turbulent  axisymmetric  flow 
patterns,  and  the  heat  transfer  to  the  tube  wall. 

Since  an  expansion  flow  is  quite  removed  from  the  phenomena  occurring 
during  a  ballistic  cycle,  a  second  type  of  flow  is  simulated  which  has  time- 
dependent  pressure  and  temperature  profiles  similar  to  a  real  weapon.  It  is 
obtained  by  adding  proper  heat  and  mass  to  the  one-phase  flow  via  source 
terms.  An  empirical  burning  law  for  pressure-dependent  sources  is  used.  Tin- 
sources  move  with  the  fLow.  We  designate  this  idealization  an  the  "real 
gun”.  The  essential  parameters  for  the  real  gun  differ  from  Table  1  mlv  wit  I 
respect  to  the  initial  conditions.  The  initial  gas  pressure  is  assumed  to  he 
ambient  pressure  (0.1  MPa)  and  the  initial  gas  temperature  to  he  ambient 
temperature  (293  K) .  In  all  cases  involving  the  real  gun  simulation,  the 
projectile  Is  released  from  Its  initial  position  when  the  pressure  at  the 
projectile  base  reaches  30  MPa. 

Because  the  DELTA  code  is  based  on  an  implicit  finite  diftorenoe  sohe;m  ■ , 
no  stability  condition  restricts  the  size  of  the  time  step.  The  pre*--,  >ted 
results  are  all  calculated  using  a  constant  time-step  of  1  1  (.s  .  The 
computational  mesh  consists  of  49  uniformly  or  nonun i formlv  spiced  mesh  point- 
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in  the  axial  direction  and  19  nonunif orraly  spaced  mesh  points  in  the  radial 
direction.  To  obtain  a  finer  spatial  resolution  in  the  boundary  region,  the 
mesh  points  are  concentrated  near  the  wall.  An  example  of  a  19x49 
computational  mesh,  which  is  used  in  most  of  the  computations,  is  shown  in 
Figure  1.  The  smallest  grid  size  in  radial  direction  at  the  bore  surface  is 
7.7  ,im  .  The  mesh  for  computing  the  heat  conduction  in  the  tube  wall  is 
generated  in  the  same  way  with  the  same  mesh  distribution  in  axial  direction 
as  on  the  gas  side,  and  a  corresponding  mesh  concentration  near  the  inner  bore 
surface.  Both  the  size  of  time-step  and  the  number  of  grid  points  seem  to  be 
reasonable  compromises  between  accuracy  and  computing  time. 

Some  of  the  computed  results  for  the  Lagrange  and  real  gun,  respectively, 
with  laminar  flows  and  adiabatic  boundaries  are  compared  in  Figures  2-5.  The 
figures  show  the  histories  of  the  gas  pressures  and  temperatures  at  the  center 
of  the  breech  and  projectile,  and  the  velocities  and  displacements  of  the 
projectile.  The  main  differences  are  in  the  temporal  distributions  of  all 
quantities  shown,  and  between  the  final  values  of  the  muzzle  velocities.  In 
the  real  gun  simulations,  the  pressure  at  the  projectile  base  reaches  30  MPa 
at  about  2.3  ms,  at  which  time  the  projectile  begins  to  accelerate  down  the 
tube . 
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figure  1.  Standard  Computational  Mesh  19  x  49.  Minimum  Ladiai 
Crid  Size  is  7.7  Urn. 


Figure  2.  Pressure  Histories  at  the  Center  of  the  Breech  and  the 
Projectile  Base  for  the  Laminar  Flow  with  the  Adiabatic 
Walls  in  the  Lagrange  Gun  (LG)  and  the  Real  Gun  (RG) 
Simulations. 
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Detailed  qualitative  resuLts  for  a  Lagrange  gun  with  a  laminar  flow  and 
adiabatic  walls  are  presented  in  Figures  6-10.  Notice  that  the  three- 
dimensional  surfaces  are  shown  from  different  directions  in  order  to  make 
their  display  clearer.  The  spatial  distribution  of  the  axial  velocity  is 
shown  in  Figures  6  and  7  at  3.75  ms.  At  this  time,  the  projectile  exits  the 
gun  tube  with  a  muzzle  velocity  of  623  m/s.  Figure  6  shows  the  axial  velocity 
field  when  a  uniformly  axially  spaced  mesh  as  in  Figure  1  is  used.  Figure  7 
shows  the  same  quantity  but  computed  using  a  nonuniformly  spaced  axial  mesh. 
The  radial  distribution  of  the  mesh  points  are  the  same  in  both  figures.  For 
a  given  axial  position  the  axial  velocity  is  constant  across  much  of  the 
radius  of  the  tube  (the  core  flow  region),  and  decreases  to  zero  only  very 
close  to  the  wall  (the  boundary  layer  region).  The  boundary  layer  is  the 
result  of  the  no-slip  condition  (w=0)  at  the  wall.  The  thickness  of  the 
boundary  layer  is  approximately  0.2  to  0.3  mm.  In  the  axial  direction,  the 
velocity  is  distributed  linearly  in  the  core  region  between  the  zero  value  at 
the  breech  and  the  muzzle  velocity.  The  three-dimensional  temperature  and 
pressure  distributions  are  given  in  Figures  8  and  9.  Due  to  the  adiabatic 
boundary  condition  at  the  wall  surface  the  heat  generated  by  the  viscous 
forces  near  the  tube  wall  cannot  transfer  to  the  tube  wall,  and  the  gas 
temperature  rises  towards  the  wall  surface.  Here  again  the  boundary  layer  is 
only  0.2  to  0.3  mm  thick.  The  pressure,  however,  stays  constant  in  radial 
direction  over  the  entire  cross-section.  The  assumption  in  boundary  layer 
theory  that  the  radial  pressure  gradient  is  zero  would  be  valid  in  this 
example.  Figure  10  shows  the  3-D  graph  of  the  radial  gas  velocity.  In 
approximately  the  first  70%  of  the  distance  to  the  projectile,  the  radial 
velocity  is  negative,  i.e.,  the  flow  is  directed  towards  the  center  line  (C- 
L) .  Thereafter,  it  is  positive  and  increases  remarkably  towards  the 
projectile  base.  Only  very  close  to  the  projectile  does  the  value  of  the 
radial  velocity  drop  rapidly  from  its  maximum  value  to  zero.  Of  course,  the 
radial  velocity  is  small  in  comparison  to  the  axial  velocity  since  it  is 
induced  only  by  molecular  viscosity  and  heat  conductivity.  The  results  in 
Figures  2-10  agree  very  well  with  both  a  one-dimensional  solution  of  the  core 
flow  using  the  method  of  characteristics,  and  the  two-dimensional  numerical 
calculations  of  Heiser  and  Hansel. 
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Figure  b.  Lagrange  Gun,  Laminar  Flow,  Adiabatic  Walls:  Spatial 
Distribution  of  the  Axial  Gas  Velocity  at  the  Time  of 
Muzzle  Clearance.  Mesh  is  Uniformly  Spaced  in  Axial 
Direction  and  Nonuni fomilv  Spaced  in  Radial  Direction 
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For  the  Lagrange  gun  considered  above  the  Reynolds  number  based  on  the 
tube's  diameter  and  muzzle  velocity  is  of  the  order  of  107.  Therefore,  the 

assumption  of  a  laminar  flow  might  not  be  realistic,  and  we  examined  in 

further  calculations  the  effects  of  turbulence.  In  comparison  to  the  laminar 
flow,  several  important  differences  were  observed  when  the  two  algebraic 
turbulence  models  described  in  Section  IV  were  applied.  However,  the 

difference  between  the  effects  of  the  two  turbulence  models  was  found  to  be 

insignif leant.  Therefore,  we  do  not  specify  the  type  of  model  in  the 
presented  results.  First,  one  observes  a  difference  in  the  projectile 
performance  between  the  laminar  and  turbulent  type  flows.  This  is  shown  in 
Figure  11  in  terras  of  the  pressure  histories  at  the  breech,  in  Figure  12  in 
terras  of  the  temperature  histories  at  the  breech,  and  in  Figure  13  in  terms  of 
the  projectile  velocity  histories.  In  the  turbulent  flow  simulation,  the 
projectile  muzzle  velocity  is  about  50  m/s  less  than  in  the  laminar  flow 
simulation.  The  axial  velocity  flow  field  at  the  time  of  muzzle  clearance  is 
shown  in  Figure  14.  The  velocity  boundary  layer  is  fully  developed  between 
the  center  line  and  the  tube  wall.  The  axial  velocity  overshoots  the 
projectile  velocity  near  the  center  line.  This  overshoot  is  related  to  the 
radial  gas  velocity  (Figure  15)  which,  near  the  projectile  base,  is  <<ne  to  two 
orders  of  magnitude  larger  than  in  the  case  of  a  laminar  flow  (Figure  10).  In 
this  region,  the  greater  radial  gas  flux  toward  the  tube  wall  transports  mass 
away  from  center-line  which  in  turn  can  accelerate  the  axial  flow.  The  radial 
variation  of  the  axial  velocity  (Figure  16),  the  radial  velocity  (Figure  17) 
and  the  temperature  (Figure  13)  taken  0.25  m  upstream  of  the  muzzle  show  some 
details  of  the  differences  between  the  laminar  and  turbulent  flows  at  the  time 
of  muzzle  clearance.  Boundary  layer  calculations  in  Ref.  7  show  comparable 
trends  between  laminar  and  turbulent  flows.  However,  the  flow  patterns 
computed  by  different  types  of  turbulence  models,  e.g. ,  non-algebraic  models, 
may  differ.  An  experiment  corresponding  to  this  idealized  expansion  flow  is 
needed  to  validate  a  turbulence  model.  Such  experiments  are  being  attempted 
at  the  French-German  Institute  (ISL)  in  France  and  the  Ernst-Mach- Institut 
(EMI)  in  Germany. 

In  a  second  series  of  calculations  we  investigated  the  significance  of 
heat  transfer  to  the  tube  wall.  The  heating  of  the  wall's  surface  has 
practical  implications  because  it  influences  the  erosion  of  the  tube.  Our 
heat  transfer  model  together  with  the  calculation  of  the  heat  conduction  in 
the  wall  couples  the  unsteady  behavior  in  both  media,  and  it  is  discussed  in 
Section  III.  We  now  show  results  obtained  for  the  laminar  Lagrange  gun 
expansion  flow  with  heat  transfer  from  the  gas  to  the  tube  wall.  Initially 
the  wall  is  assumed  to  be  at  ambient  temperature.  The  thermal  properties  of 
the  gun  barrel  are  characterized  by  the  barrel  material  density  p  =7.8  kg/m  , 
thermal  conductivity  \  =43  W/(m  K)  and  specific  heat  cw  =  460  J7(kg  K).  The 

pressure  and  temperature  histories  at  the  breech,  and  projectile  velocity 
history  are  given  in  Figures  11-13,  and  can  be  compared  with  the  other  two 
simulations.  The  spatial  profile  of  the  axial  velocity,  radial  velocity, 
pressure  and  temperature  are  plotted  in  Figures  19-22,  respectively,  at  the 
time  of  muzzle  clearance.  Comparing  the  results  to  those  of  the  laminar  flow 
with  adiabatic  walls,  we  find  that  the  velocity  boundary  layer  as  well  as  the 
temperature  boundary  layer  are  in  the  present  case  thinner,  that  the  radial 


Figure  12.  Lagrange  Gun:  Temperature  Histories  at  the  Center 
of  the  Breech  for  Both  Laminar  and  •Turbulent  Flows 
with  Adiabatic  Walls,  and  for  Laminar  Flow  with 
Heat  Transfer. 
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Figure  13.  Lagrange  Cun:  Projectile  Velocity  Histories  for  Loth 
Laminar  and  Turbulent  Flows  with  Adiabatic  Walls  and 
for  Laminar  Flow  with  Transfer. 
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Figure  16.  Lagrange  Gun,  Adiabatic  Walls:  Radial  Profiles  of  the 
Axial  Gas  Velocity  for  Both  Laminar  and  Turbulent  Flow 
at  the  Time  of  Muzzle  Clearance  at  0.23  m  Away  from 
the  Muzzle. 
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Figure  17. 


Lagrange  Cun,  Adiabatic  Walls:  Radial  Profiles  of 
Radial  Cas  Velocity  for  Both  Laminar  and  Turbulent 
Flow  at  the  Time  of  Muzzle  Clearance  at  1.25  m  Away 
from  the  Muzzle. 
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Figure  18.  Lagrange  Gun,  Adiabatic  Walls:  Radial  Profiles  c' 
the  Cas  Temperature  fcr  Both  Laminar  and  1  urbu  1  out. 
Flows  at  the  Time-  of  Muzzle  Clearance  at  C  . 25  n- 
Away  from  the  Muzzle. 


Lagrange  Gun,  Laminar  Flow,  Heat  Transfer:  Spatial 
Profile  of  the  Axial  Cas  Velocity  at  the  lime  of 
Muzzle  Clearance. 
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velocity  is  about  one  order  of  magnitude  higher,  and  that  the  radial  velocity 
is  always  directed  toward  the  wall.  For  example,  the  velocity  boundary  layer 
is  only  about  0.05  mm  thick.  Detailed  comparisons  across  the  boundary  layer 
regions  for  axial  velocities  and  temperatures  are  presented  in  Figures  23  and 
24,  and  across  the  tube  cross-section  for  the  radial  velocities  are  shown  in 
Figure  25.  All  these  profiles  are  at  0.25  m  from  the  muzzle  at  the  time  of 
muzzle  clearance.  Figure  26  shows  the  history  of  the  wall  surface  temperature 
at  two  different  locations  along  the  tube  wall.  One  is  taken  at  150  mm  away 
from  the  breech.  This  wall  point  belongs  to  the  chamber,  and  is  heated  up 
from  the  beginning  of  the  flow  cycle.  The  second  wall  point  is  250  mm  away 
from  the  breech,  and  is  heated  up  only  after  the  projectile  has  passed  it  (at 
0.47  ms).  The  maximum  wall  temperature  occurs  early  in  the  4  ms  cycle.  The 
laminar  gas  layer  close  to  the  tube  wall  cools  rapidly  because  the  transport 
of  heat  in  the  tube  wall  by  conduction  is  much  faster  than  the  transport  of 
heat  by  conduction  and  convection  on  the  gas  side  towards  the  wall.  To  obtain 
the  heat  transfer  precisely  from  the  gas  to  the  wall,  we  need  to  compute  the 
radial  temperature  gradient  on  both  sides  of  the  inner  bore  surface  as 
accurately  as  possible.  The  temperature  boundary  layer,  however,  is  very  thin 
as  it  is  shown  in  Figures  22  and  24.  This  implies  that  very  small  grid  sizes 
must  be  used  close  to  the  wall.  The  smallest  one  in  radial  direction  for 
Figure  26  is  0.73  pm  . 

The  final  simulations  are  of  one-phase  flows  in  the  real  gun,  i.e.,  of  a 
gas  flow  with  mass  and  energy  sources.  The  pressure  and  temperature  histories 
at  the  breech  and  projectile,  and  the  projectile  velocity  and  displacement 
histories  for  the  laminar  flow  simulation  and  for  the  real  gun  are  given  in 
Figures  2-5.  Figures  27-40  are  a  series  of  three-dimensional  profiles  which 
characterize  both  the  laminar  and  turbulent  flows  inside  an  adiabatic  tube. 

The  group  of  figures  associated  with  the  laminar  flow  are  at  two  times:  the 
first  at  3.6  ms  when  maximum  pressure  is  achieved  and  the  second  at  5.3  ms 
when  the  projectile  exits  the  cube.  For  these  calculations  a  nonuniforra 
radial  grid  is  used  with  the  minimum  grid  size  of  0.73  ym  because  of  the 
thinness  of  the  boundary  layer.  We  excluded  the  profile  of  the  radial 
velocity  at  3.6  ms  because  the  magnitude  of  this  component  was  smaller  than 
0.05  m/s.  Three-dimensional  profiles  show  the  turbulent  flow  fields  when 
maximum  pressure  occurs  at  3.6  ms  (Figures  34-37),  and  when  the  projectile 
exits  the  tube  at  5.49  ms  (Figures  38-40).  As  in  the  Lagrange  gun,  the 
differences  are  minimal  between  the  simulations  with  different  turbulence 
model,  but  significant  between  the  laminar  and  turbulence  simulations. 
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Figure  24.  Lagrange  Gun,  Laminar  Flow:  Radial  Profiles  of  the 

Gas  Temperature  for  Adiabatic  and  Heat  Termeable  Wall: 
at  the  Time  of  Muzzle  Clearance  0.25  m  Away  from  the 
Muzzle. 
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Figure  36. 


Real  Gun,  Turbulent  1 1  ov, ,  Adiabatic  Wallet 
Spatial  Distribution  of  the  Gas  Temperature 
at  3.6  ms. 
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Hgure  W.  Real  Gun,  Turbulent  Flow,  Adiabatic  Walls: 
Spatial  Distribution  of  the  Axial  Gas 
Velocity  at  Muzzle  Clearance  (5.49  ms). 


igure  40.  Real  Gun,  Turbulent  Flov:,  Adiabatic  Walls: 

Spatial  Distribution  of  the  Gas  Pressure  at 
Muzzle  Clearance  (5.49  ms). 


VI.  SUMMARY 


Special  interior  ballistic  phenomena  are  investigated  by  using  the  DELTA 
computer  code,  which  is  designed  to  calculate  the  multidimensional,  two-phase 
flow  behind  an  accelerating  projectile  inside  a  gun  tube.  Details  of  the  heat 
transfer  and  turbulence  submodels  used  in  the  simulations  are  given  in  this 
report.  The  general  mathematical  model  and  numerical  algorithm  are  described 
in  a  companion  paper,  Ref.  3.  Results  are  given  for  two  types  of  idealized, 
one-phase  interior  ballistics  gun  simulations:  a  pure  expansion  flow,  and  a 
flow  with  moving  mass  and  heat  sources.  Comparisons  are  made  between  laminar 
and  turbulent  flows  as  well  as  between  flows  in  an  adiabatic  tube  and  in  a 
tube  that  allows  heat  transfer.  The  comparisons  show  that  the  flow  structure 
and  the  motion  of  the  projectile  are  both  significantly  influenced  by 
turbulence  and  heat  conduction,  because  of  a  strong  coupling  between  the  core 
flow  and  boundary  phenomena.  A  study  of  the  effects  of  turbulence  and  heat 
conduction  therefore  can  be  best  accomplished  with  a  multi-dimensional 
mathematical  model,  such  as  the  DELTA  code. 
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